APS(Hadron)/AYNU[20 13] 
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An isoscalar resonance, X(1870), with a mass M = 1877.3 ± 63{stat)*j'l{syst)MeY and width F = 57 ± 
12(.sfa/)+^'(i'y.sf)MeV was observed by the BesIII Collaboration recently. The mass of X( 1870) is compatible 
with 1/2(1 870) but the decay width is much smaller. In this paper, we analyzed the mass and the two-body 
strong decay of 772(1870) and X(1870). According to the mass spectra predicted by the Godfrey-Isgur potential 
model, ^(1870) can be candidates of the states f[{2^Pi), ri'^il^Di) or riiil^Di). In the Regge trajectories (RTs) 
framework, X(1870) favors the ;7'(3'5o) and f[{'2?Pi) assignments. However, the decay properties of ;7'(3'5o), 
f[{2^P\), ;/;i(l'£)2) disagree with the experiments. The two-body strong decay widths indicate that 7/2(1870) and 
X(1870) maybe the two different resonant structures. 7/2(1870) could be assigned to the r}2{2^D2) quarkonium 
as the first radical excited state of 7/2(1645). X(1870) favors a 7/;/(0^2"^) hybrid assignment. Some significant 
decay modes are suggested to test in future. 

PACS numbers: 12.38.Lg, 13.25.Jx 



I. INTRODUCTION 

Quantum Chromodynamics (QCD) not only allows the qq 
mesons to exist but also enables the exotic states, e.g., glue- 
balls, hybrids, multi-quark states, molecules, ..., to appear in 
the nature. Although people paid a lot of efforts, no unique 
experimental signatures of such exotic states have been well 
established yet jTI]. 

For hybrids, they can carry exotic J^'- quantum numbers 
which are not accessible by the conventional qq mesons in 
the constitute quark model. Experimentally, there are several 
candidates with exotic J^^ ^. Among them, 7ri(1400) and 
7ri(1600) maybe the most famous states. Due to its low mass 
and a single decay mode, ;ri(1400) looks more like a multi- 
quark state rather than a hybrid meson. Identification of the 
7ri(1600) as the lightest hybrid seems also premature because 
the current measurements are not inconsistent with a multi- 
quark interpretation [2J . 

Hunting the hybrid mesons with the canonical 7^'' seems 
more difficult for the following reasons. Firstly, the predicted 
mass spectra of hybrid in different models usually overlap the 
high-excited states of quarkonium. Secondly, both excited 
quarkonium and hybrid can decay into L - and L - \ 
("5 -I- P") qq meson pairs. Finally, the mixing will occur when 
they carry the same y^*-. 

A 1 strangeonium-like resonance, 7(2175), which was 
observed by BaBar Collaboration |3, 4] and confirmed by the 
BES |5] and Belle la] Collaborations, was interpreted as a 
ssg candidate [31 ■ However, the other assignments are still 
possible, e.g., a 2^D\ ss state (a], a ssss tetraquark state ^ 
[Toll , or a threshold effect [11]. The other hybrid candidates, 
e.g., the 2"^ states of ;r2(1880) and 772(1870), have never been 
established well. 



It is particularly interesting that 772(1870) has been observed 
in yy reacti on IIT2I fl3ll . radiative Jjip decays Ill4ll and pp 
annihilation jlSlflql . Its mass determined in these processes 
are about 1 870]VIeV and the decay width was confirmed to be 
large (around 200]VIeV). Now its available mass and width 
listed by Particle Data Group (PDG) are 



M = 1842 ± 8]VIeV, 



r = 225 ± 14MeV. 



It should be stressed that the radiative y/i/r decay ( FiglTtA]) 
and pp annihilation are the ideal glueball hunting grounds. 
But glueball production is suppressed in yy reaction. 



j/iii 
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mmxxm 
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FIG. 1: [A]. A prior production of glueball in the J lip — > -yXc; [B]. 
A prior production of hybrid in the Jjij/ — > wXh- 



A new resonant structure of X(1870) was observed by 
BesIII with a statistical significance of 7.2cr in the channels 
7/i/f -^ (jj + X(1870) — > cijrin^TT^ recently 111711 . Its mass and 
width were given as 



M= 1877.3: 



: 6.3+^;4]VIeV, 



r = 57 + 



12+^'']VIeV. 
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Here the first errors are statistical and the second ones 
are systematic. The mass of the X(1870) is consistent with 
7/2(1870) but the width is much smaller. By contrast, X(1870) 
was observed in the hadronic J/tfr decay which is considered 
"hybridrich"(FiginB]). If the process, J/tfr -^ w-i-X(1870),is 



mainly hadronic, the quantum numbers of X( 1870) should be 
0+0"+, 0+1++ or 0+2-+. The predicted masses for the 0+0"+, 
0+1 ++ or 0+2+ hybrids overlap 1.8GeV(see in Table HV) in 
different models, e.g., the Bag model il8l[l9ll . the flux-tube 
model [20., 211 . the constituent gluon model ||22I1 . Hence, 
X(1870) becomes a good hybrid candidate. 

At present, 772(1870) with 7^*- = 2"+ is listed as a member 
of 1 'D2 nonet in PDG with 772(1645), 712(1670) and KiimO)' 
temporarily [1]. 772(1 870) was only observed in the flo(980)7r, 
72(1270)77 , and fl2(1320)7r channels, not in the KKn channel, 
then it is hard to regard 772(1870) as an excited ri'(ss) state 
unless large nil <-> ss {nh = {uu + dd)/2) flavor mixing in the 
772(1870). However, the near degeneracy of the ?72(1645) and 
;7r2(1670) argue against this mixing possibility. 

772(1870) has been examined as a quarkonium in the ■'Po 
model [23] . The authors argued that the interpretation of 
772(1870) as a 2^D2(nn) state is possible. However, the 
hybrid interpretation was not excluded there. In addition, 
the flo(980)7r channel which has been observed for 772(1870) 
wasn't calculated in Ref. [23]. 

So it is necessary to exhaust the nature of 772(1870) and 
X(1870) systematically. Some crucial tests to distinguish 
them are also necessary. In this paper, we will study the 
masses and strong decays of X( 1870) and 772(1870) in detail. 
The paper is organized as follows. In Sec. II , the masses of 
X(1870) and 772(1870) are explored in the relativized quark 
model and the the RTs framework. In Sec. Ill , the decay prop- 
erties of these excited isoscalar quarkonium will discussed 
where the SU(3) favor symmetry has been assumed. Then 
the two-body strong decay of these states will be explored 
within the ^Pq model and the flux-tube model. Finally, we 
will present our conclusions and discussions in Sec. IV. 



II. MASS SPECTRUM 



leading-order perturbations which include the spin-spin con- 
tact hyperfine interaction, spin-orbit and tensor interactions 
and a longer-ranged inverted spin-orbit term. They arise from 
one gluon exchange (OGE) forces and the assumed Lorentz 
scalar confinement. The expressions for these terms may be 
found in Ref H. 

Now, the spin-averaged mass, M„i, of riL multiplet can be 
obtained by solving the spinless Salpeter equation which is an 
approximation of the Bethe-Salpeter equation 11261 1271] . 



[^pl+ml + 



i 



Pl + ml + V{r)]ilj{r) = Eijjir). (3) 



Here we employ a variational approach described in 
Ref. Vm to solve the Eq.(3). This variational approach has 
been applied in solving the Salpeter equation for cs ||29|] . cc 
and bb [30] mass spectrum. 

In the calculations, the simple harmonic oscillator (SHO) 
functions are taken as the trial wave functions. It is given by 



^„t(r,/3)^/3- 



3/2 



2(27? - 1)! 
'r(« + /+i) 



'+l/2/'«2 2x 



Oer)'e-— L;iY^(^^rO 



in the position space. Here the SHO function scale /3 is the 
variational parameter 

By the Fourier transform, the basic SHO radial wave func- 
tion in the momentum is 



(-1)" 2(2n-l)\p, ^^ .,,i.p\ 



0^1^ \Y{n + I + W P' 



/3' 



The wave functions of if/„i(r,/S) and 4'ni(p,l3) meet the nor- 
malization conditions: 

iA>,/?)r2fl'r=l; il,l,(p,p)p^dp ^ 1. 

Jo 



Since the masses of X( 1870) and 772(1870) are nearly equal, 
only X(1870) will be discussed in this section for convenience. 
Obviously, the conclusions are also suitable for 772(1870). 

In the Godfrey-Isgur relativized potential model [24], the 
Hamiltonian consists of the central potential and a kinetic term 
in a "relativized" form 



In the variational approach, the corresponding M„i are given 
by minimizing the expectation value of H 



where 



dp 



(4) 



H- ^Pl, + ml+ ^p\ + ml + V{r). (1) 

The funnel-shaped potentials which includes a color 
coulomb term at short distances and a linear scalar confin- 
ing term at large distances are usually incorporated as the 
zeroth-order potential. The typical funnel-shaped potential 
was proposed by the Cornell group (Cornell potential) with 
the form |125D 



4 a 

yqqir) = --— +o-r + C 

i r 



(2) 



The strong coupling constant o-j, the string tension a- and 
the constant C are the model parameters which can be fixed 
by the well established experimental states. The remaining 
spin-dependent terms for mass shifts are usually treated as the 



Enlifi) = {H)„, = {ljJ„l\H\iP„i). 



(5) 



When all the parameters of the potential model are known, 
the values of the harmonic oscillator parameter^ can be fixed 
directly. With the values of y6, all the spin-averaged mass M„i 
will be obtained easily. M„i obtained in this way trend to be 
better for the higher-excited states [31]. 

It is unreasonable to treat the spin-spin contact hyperfine 
interaction as a perturbation for the ground states, because 
the mass splitting between pseudoscalar mesons and vector 
mesons are much large. Then we consider the contributions 
of ^sK'") for the 15 mesons. The following Gaussian-smeared 
contact hyperfine interaction 0321] is taken for convenience. 



9'«o ^171 



(6) 



In this work, we choose the model parameters as follows: 
OT„ = OTrf = 0.220 GeV, m, = 0.428 GeV, c, = 0.6, a = 0.143 
GeV^, K = 0.37 GeV, and C = -0.37 GeV. A smaller cr was 
obtained by the relation between the slope of the Regge trajec- 
tory for the Salpeter equation a' and the slope a'^, in the string 
picture 03311 . This smaller value for cr is taken here rather than 
the value in Ref. \M\. The Gaussian smearing parameter k 
seems a little smaller than that in Ref. |24]. However, the k 
is usually fitted by the hyperfine splitting of low-excited nS 
states in the literatures with a certain arbitrariness. 

The values of M„l and p for the states 25 , 35 , 45 , IP, 2P, 
3P, ID, 2D, 3D, IF, IF, IG and 1// ai-e listed in Tableffl The 
experimental masses for the relative mesons are taken from 
PDG|1]. 



States 


M„i{nn) 


P 


M,„UJ 


M„,(i.s) 


P 


M„,UJ 


IS 


- 


0.44 0.34 


- 


- 


0.42 0.39 


- 


2S 


1.399 


0.310 


1.389 


1.631 


0.330 


1.629 


3S 


1.859 


0.295 




2.069 


0.310 




4S 


2.240 


0.290 




2.436 


0.300 




IP 


1.252 


0.310 


1.257 


1.460 


0.340 


1.478 


2P 


1.711 


0.294 




1.926 


0.315 




3P 


2.110 


0.290 




2.308 


0.300 




ID 


1.661 


0.280 


1.672 


1.883 


0.300 




2D 


2.067 


0.276 




2.272 


0.292 




3D 


2.417 


0.275 




2.609 


0.288 




IF 


1.924 


0.277 




2.128 


0.295 




2F 


2.287 


0.275 




2.478 


0.290 




IG 


2.161 


0.275 




2.350 


0.292 




IH 


2.377 


0.273 




2.554 


0.287 





TABLE I: The spin-averaged mass (unit: GeV) and the harmonic 
oscillator parameter p (unit: GeV"') of the states 25, 35, 45, IP, 
IP, 3P, ID, 2D, 3D, IF, 2F, IG and IH. 

Obviously, the spin-averaged masses of the 25, IP, ID nh 
and IP, 25 ss mesons are consistent with the experimental 
data. Indeed, the predicted masses of higher excited states 
here are also reasonable, e.g., 04(2040) and /4(2050) are very 
possible the P-wave nh isovector and isoscalar mesons with 
the mass 1996!,^"MeV,2018±llMeV Q]. The predicted spin- 
averaged mass of IP is compatible with experiments. The 
presented M„l in Table U are also overall in good agreement 
with the expectations from Ref. |34]. 

In the calculations, we ignored the Coulomb interaction for 
ID, 2D, IP, IG and IH states. In this way, M„i for these states 
increase about lOOMeV As shown in Ref. 134,13511 . the masses 
for higher excited states with cr = 0.143GeV^ and Os - are 
closer to experimental data than the results given in Ref. 02411 . 

The trend that a higher excited state corresponds to a 
smaller y6 coincides with Ref. |36-38]. It needs to mention 
that the strong decay amplitudes given by the ^Pq model often 
contain the nodal-type Gaussian form factors which can lead 
to a dynamic suppression for some channels. Then the values 
of j6 are important when one use the ^Pq model to exact the 
decay width for the higher excited mesons. 



The masses of 77'(3>5o), /i'(2^Pi), t/^CI'Dz) and ri2{2^D2) 
are usually within 1.9 ~ 2.1 GeV in various quark potential 
models (see in Table Hlb. The predicted regions are compati- 
ble with the spin-averaged masses of 3S{ss), 2P{ss), \D{ss) 
and 2D(nn) presented in Table |I] Due to the uncertainty of the 
potential models, absolute deviation from experimental data 
are usually about 100~150MeV for the higher excited states. 
Comparing with these predicted masses, X(1870) disfavor the 
77'(3'5()) assignment for its low mass. But the possibilities of 
f[{2?Pi), /72(1'D2) and 772(2'D2) still exist. 



States 



Ref. (24|] 
Ref. m 
Ref. liSl 



Y(3'5o) /;(2V,) 7/^(1 '£>2) m&D.) 



2030 
2085 2016 

2099 1988 



1890 2130'' 

1909 1960 

1851 



TABLE II: The masses predicted for 3'5o(77'), 2^,(77'), l^Diirj') and 
2'Z)2(?7)inRef. |24, 39, 40]. 

Regge trajectories (RTs) is another useful tool for studying 
the mass spectrum of mesons. The linearity and parallel of 
RTs in heavy-light and heavy-heavy mesons is still in debate. 
In Ref. 14 111 , the authors fitted the RTs for all light-quark me- 
son states listed in the PDG tables. A global description was 
constructed as 



M^ = 1.38(4)n+ 1.12(4)7- 1.25(4). 



(7) 



/zi(1380), /i(1420) and ?7'(1475) have been established as 
the I'Pi, 1^1 and 2'5o ss mesons in PDG |1]. With the 
mass difference between the mass squared of X(1870) and 
these states (Table HIH). X(1870) could be assigned for the 
77'(3'5o) and f[{2^Px). The mass of X(1870) is too large for 
the 772(1 'D2) state in the RTs. 7/2(1645) has been assigned as 
the 1 ID2 nh meson fl. Since m2(X(1870)) - m2(?72(1640)) * 
0.91 which is much smaller than 1.38(4), X(1870) can't be 
accepted as the 2'D2 tjh state for its low mass. 

The mass of X( 1870) in the RTs 



77'(3'5o) 
77'(1475) 
1.34 



nil' Pi) 
/i(1420) 
1.38 



77^(1'D,) 
/7i(1380) 
1.60 



772(2'£)2) 

772(1645) 
0.91 



TABLE III: X(1870) calculated in RTs for different states are shown. 

As mentioned in the Introduction, X(1870) is also a good 
hybrid candidate since its mass overlaps the predictions given 
by different models (Table ITVll. 



States 

Bag |i8|,lig|] 



Flux tube |20|,|21D 
Constituent gluon (22 



77^,(0+0-+) /«(0+r+) 77«(0+2-+) 



1.3 

1.7-1.9 

1.8-2.2 



heavier 
1.7-1.9 
1.3-1.8 



1.9 

1.7-1.9 

1.8-2.2 



TABLE IV: The masses predicted for 77^(0+0"+) ///(O+r 
7/«(0+2-+) hybrid states in Ref. L18r,22i. 



) and 



In summary, the X(1870) could be a f^il^Pi), r]'^(l^D2) or 
772(2 'D2) meson in the Godfrey-Isgur quark potential model 
with a reasonable uncertainty. In the RTs, X(1870) favors the 
?7(3'5'o) and f'^{2^P\) assignments. X(1870) is also a good 
hybrid state candidate. The investigations of the strong decay 
properties will be more helpful to shed light on the underlying 
structure of this state. 



III. THE STRONG DECAY 
A. The final mesons include the scalar mesons below 1 GeV 

Despite many theoretical efforts, the scalar nonet of light 
mesons has never well-established. The lowest-lying scalar 
mesons including cr(500) (or /o(600)), a-(800), flo(980) and 
/o(980) are difficult to be described as qq states, e.g., flo(980) 
is associated with nonstrange quarks in the qq scheme. If this 
is true, its high mass and decay properties are difficult to be 
understood simultaneously. So interpretations as exotic states 
were triggered, i.e., as two clusters of two quarks and two 
antiquarks [42], particular quasimolecular states [43], andun- 
correlated four quark states qqqq [44] have been proposed. 

Though the structures of these scalar mesons below 1 GeV 
are still in dispute, the viewpoint that these scalar mesons can 
constitute a complete nonet states has been reached in most 
literatures (as illustrated in Fig. 2). In the following, we will 
denote this nonet as " S " multiplet for convenience. 
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FIG. 2: The S nenet below 1 GeV shown mY - h plane. 



Due to the unclear mechanism, it seems much difficult to 
study the decay processes where the final mesons includes a S 
member. In Ref. [45, 46], ao(980), cr(500) and /o(980) were 
treated as I^Pq qq mesons approximately. In Ref 11231 14711 . 
the decay channels including the S meson have been ignored. 
However, this kind of decay channel maybe predominant for 
some mesons. For example, the observations indicate that 
X(1870) primarily decay via the flo(980)7r channel [ 17]. 

In what follows we will obtain some useful information 
about this kind of decay mode where the SU(3) flavor symme- 
try shall be invoked. We will show that flo(980)7r, o-qt] and /qt; 



are the main decay channels for the isoscalar nh and the nhg 
mesons when they decay primarily through " ^S + P" mesons, 
where the sign "P" denotes a light pseudoscalar meson. This 
explains why X(1870) has been first observed in the rjn^n^ 
channel. 

We noticed that the S nonet could be interpreted Uke the 
qq nonet in the diquark-antidiquark scenario. In Wilczek and 
Jaffe's terminology [48, 49], the S mesons consist of a "good" 
diquark and a "good" antidiquark. When 11, d quarks forms 
a "good" diquark, it means that the two light quarks could 
be treated as a quasiparticle in color 3, flavor 3 and the spin 
singlet, which denotes as [ud\. 

In the diquark-antidiquark lirnitj the parity of a tetraquark 
is determined by P = (-1)^'-''^ 15(| where the L12-34 refer to 
the relative angular momentum between two clusters. Thus 
the S mesons are the lightest tetraquark states in the diquark- 
antidiquark model with L12-34 - 0. The S nonet in the full set 
of flavor representations is 

(3 (g) 3)3 (g) (3 ® 3)3 = 8el 

Because the SU(3) flavor symmetry is not exact, two phys- 
ical isoscalar mesons, cr() and^ are usually the mixing states 
of the |8)/=o and |l)/=o states 114211 . 



/o 

CTO 



COS)? sini?W|8>/=o 
-sini? COS)? )\ |l)/=o 



(8) 



When the mixing angle d- equals the so-called ideal mixing 
angle, i.e., 1? - 54.74°, the composition of the cr(500) and 
/o(980) are 



/o 
a-{) 



_ I \-^{[su\[su\ 



sd\[sd\)) 
\[ud\[ud'\) 



The deviation from the ideal mixing angle of the the cr(500) 
and /o(980) is small 14211 . In the following calculations, we 
will treat them in the ideal mixing scheme. 

Under the SU(3) flavor assumption, all the members of the 
octet have the same basic coupling constant in one type of 
reaction, while the singlet member have a different coupling 
constant. Particularly, when a quarkonium decays into S and 
qq mesons, there are five independent coupling constants, i.e., 
gAii gAi\, gA\i, gB&a and gsn, corresponding to five diff'erent 
channels 



1 8V7 ^1 8)5® 1 8%g : 


gA&S 


1 8),^ ^1 8)^® 1 1>,^ : 


gASl 


1 8),^ ^1 1)5® 1 8V7 : 


gAlS 


1 i)gg ^1 8)5® 1 8),^ : 


gBSS 


1 l>w ^1 1>5® 1 l)w : 


gB\l 



In order to determine the relations between these coupling 
constants, we shall assume the process that a qq meson decays 
into a S and another qq mesons obeys the OZI (Okubo-Zweig- 
lizuka) rule, i.e., the two quarks in the mother meson go into 
two daughter mesons respectively. Therefore, there are four 
forbidden processes, X(-j={uli - dd)) -^> ao + ss, X{-j=(uii + 



dd)) -f^ (Tq+ss, X{ss) -f> fo+-j={uu+dd) andX{ss) ->* ctq+ss. 
With the help of the SU(3) Clebsch-Gordon coefficients IH, 

I 



the ratios between the five coupling constants are extracted by 
the four forbidden processes as 



31 : gAis. ■ gBSi» ■ gBU ■ gA&si ~ 1.41 : -2.05 : -2.89 : -2.05 : 1 

I 



(9) 
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FIG. 3: The coefficients ^^ of/ versus the mixing angle 6. 




) (Degree) 



FIG. 4: The coefficients ^^ of /' versus the mixing angle 6. 



It is well known that the physical states, 7/(548) and 77'(958) 
are the mixture of the SU(3) flavor octet and singlet. They are 
written in terms of a mixing angle, Op, as follows 



' '7(548) \ _ I 
, '7'(958) j - \ 



cos 6p - sin 9p W |8)/=o \ .-.r.. 

sinOp cos0p j\\l}i=QJ ^^' 



The mixing angle Op has been measured by various means. 
However, there is still uncertainty for Op. An excellent fit 
to the tensor meson decay widths was performed under the 
SU(3) symmetry, and Op - -17" was obtained [52]. In our 
calculation, we take Op as -17°. The excitation of nn and ss 
are denoted as 



In this scheme, the ideal mixing occurs for the choice of 
= 35.3°. When / and /' decay into a S and pseudoscalar 
mesons, the relations of decay amplitudes are governed by 
the coefficients (^ which are model-independent. With the 
coupling constants in hand, the coefficients (^ of / and /' 
versus the mixing angle are shown in the Fig |3] and Fig|4] 
When / and /' occurs in the ideal mixing, the values of (^ are 
presented in Table IV] where the mixing of 7/(548) and 7/'(958) 
has been considered. 



Decay 
channels 


aon 


err] 


kK 


fol 


O-T]' 


(-[ss] 


2.17 
0.00 


3.56 
0.00 


0.47 
0.72 


1.07 
1.08 


0.44 
0.00 



f(nn) 



cosO sinO \ / |l)/=o 
-sine cos OJX |8)/=o 



(11) 



TABLE V: .The coefficients (^ of / and /' in the ideal mixing. 

It is sure that the (^ are zero for the processes, /' — > a^TT, 
f — > err] and /' -^ ctt]', since they are OZI-forbidden. (^ of 
/' ^ /()'?' hasn't been considered in TablelVJbecause X(1870) 
lies below the threshold of /qt/'. 

As illustrated in the Fig|3] and FigH] the primary decay 
channels for f'{ss) are /ot/ and kK when the mixing angle of 
f{nn) and /'(ss) is not large. Indeed, only the ground 0"^ and 
the 0^"^ isoscalar mesons are not close to ideally mixed. Since 
X(1870) primarily decay via the ao(980);7r channel, it should 
be a (uii + dd)/ V2 predominant meson. 

What's more, 7/2(1870) was observed in the flo(980)7r, 
/2(1270)7/ , and fl2(1320)7r channels, not in the KKn chan- 
nel. Therefore, the only possible assignment for 7/2(1870) and 
X(1870) as a canonical quarkonium is the 7/2(2*^2) state. If 
7/2(1870) occupies this state, X(1870) shall be a hybrid meson 
for its narrow width [53]. Next, we will explore the two-body 
strong decay of them to pin down this possibiUty. 



B. The strong decays of 772 ( 1 870) and J?( 1 870) 

Firstly, we employ the quark-pair creation model, also 
known as the ^P^ model ll54ll55ll . to compute the two-body 
strong decay of X(1870) as the 7/2(2*^2) qq meson. This 
model has achieved great success in calculating the two-body 
strong decay of different hadrons (see Ref . ll56ll and references 
therein). 

In the nonrelativistic limit, the transition operator T of the 
^Po model is depicted as 



T = 



-lyYjX^m; 1, -m|0, 0) ff d%d%S\k3 + fe4)J/'"( 



ki - ka 



(3,4)^ (3,4), (3,4) ,t 






d^.(k3)d^ 



4; 



(^4) 



(12) 



r 



Here, co^ ' and (^J, ' are the color and flavor wave func- 
tions of the g'3^4 pair created from vacuum. Thus, Wq - 

(RR + GG + BB)I V3, ipf^^ = {uu + dd + ss)/ V3 are color and 
flavor singlets. The pair is also assumed to carry the quantum 
numbers of 0^^, suggesting that they are in a ^Pq state. Then 
x\ '_Ij represents the pair production in a spin triplet state. The 

solid harmonic polynomial }/"'{k) = \k\^'"(0ic, (pk) reflects the 
momentum-space distribution of the q^q^. 

y is a dimensionless constant which denotes the strength of 
the quark-antiquark pair created from vacuum. It is usually 
extracted by fitting the measured partial decay widths. More 
details of y will be discussed later The helicity amplitude 
j^Mj^,Mj,,Mjc(p) of A ^ B + C is given by 



{BC\t\A) = 6\Pa -Pb- Pc)M^''' ''^''> -^'^ (p), ( 1 3) 

where p represents the momentum of the outgoing meson in 
the rest frame of the meson A. When the mock state 115711 is 
adopted to describe the spatial wave function of a meson, the 
helicity amplitude M^J'^-'^^B'^^c[p) can be constructed in the 
L- S basis easily IdSII . The mock state for A meson is 



\A(n/'^-'L'f^^(P,)) 



Mtf, Ms_^ 



f dpAtp^„f'Ukuk2)\qi(k, 



)q2{h))- 



(14) 

To obtain the analytical amplitudes, the SHO wave func- 
tions are usually employed for (/'„'J '" {k\ , ^2)- For comparison 
with experiments, one obtains the partial decay width M''^{p) 
via the Jacob-Wick formula IdSII 



V2L+ 1 v^ 



2Ja 



Mj^M 



Jb,'»Jc 



(15) 



x{Jb,Mj,Jc,Mj,\JMj,}M''''^ 



.Mj„,Mjr 



(p). 



Finally, the decay width r(A -^ EC) is derived analytically 
in terms of the partial wave amplitudes 



r(A ^ BC) = 2^^^P Yj \Mls{p)\ 



(16) 



More technical details of the ^Pq model can be found in 
Ref. 1I55I1 . The inherent uncertainties of the ^Pq decay model 
itself have been discussed in the Ref. 



The nature of y has been investigated in several works. 
Firstly, the relations of the ^Pq model to "microscopic" QCD 
decay mechanisms have been studied in Ref |62]. There, the 
authors found that the constant y corresponds approximately 
to the dimensionless combination, crlmqfi, where/? means the 
meson wave function scale, niq is the mass of produced quark 
and cr is the string tension. Secondly, the momentum depen- 
dent manner of y has been studied in Ref |61]. It was found 
that y is dependent on the relative momentum of the created qq 
pair, and the form of y(k) - A + B exp(-CA:^) with k - 1^3 - £t| 
was suggested. Thirdly, J.Segovia, et al., proposed that y is 
a function of the reduced mass of quark-antiquark pair of the 
decaying meson ||63I1 . 

Based on the first and third points above, y will depend on 
the flavors of both the decaying meson and produced pairs. 
In the following, we treat the y as a free parameter and fix it 
by the 8 well-measured partial decay widths which are listed 
in TabldVIl The Mls amplitudes of these decay channels are 
presented explicitly in the Appendix A. 
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2.3 


3.91 


Pi —> jm 
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4.2 


- 


fi -^ KK 


0.579 
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5.66 



TABLE VI: . Values of 7 in different channels and comparison with 
the results given by Ref.161]. p(770), 02(1320), /;(1420), /2(1270), 
/2'(1525) andp3(1690) have been studied. 



These 8 decay channels are divided into two groups: one 
is nn -^ ns + ns, the other includes ss -^ ns + sh and 
nh —^ ns + sh. The values of 7 here are a little different from 
these given in Ref. 16 ill where an ALl potential was selected 
to determine the meson wave functions. Of course, the meson 
wave function given by different potentials will influence the 
values of y. 

It is clear in Table IVll that y decrease when/? increase. In 
addition, our calculation indicates that y depend on flavors of 
both the decaying meson and the produced quark pairs. For 
example, values of y fixed by 122 — > KK and /j — > KK are 
roughly equal. 

In the following calculations, we assume that the values of 
y corresponding to the processes of ss -^ ns + sh and nh — > 
ns + sh are determined by one function. Similarly, we take 
the function, y{p) = A+ B exp(-Cp^), for the creation vertex. 
With the four decay channels listed in fifth column of Table 
IVII we fix the function as y(p) =1.8-1- 4exp(-10/?^). For the 
processes of nh ^ nh + nh (the first column of Table IVIb. we 



fix the creation vertex function as yip) = 3.0 + 25 exp(-4;?^). 
The dependence of y on the momentum p are plotted in the 
Fig. |5] Obviously the functions can describe the dependence 
of y and p well. The functions of creation vertex given here 
need to be tested further 



n: u, d quark 



m nn^nn+nn 



nn^ns+sn 
X s~s ^ n~s + s~n 



*>., 



■■■■»- 



p (GeV) 

FIG. 5: The functions of y(p) = A + 6exp(-Cp^) in different decay 
processes. The symbols of red "•" and black "■" denote y values 
determined by the experimental data. 



Since we neglected the mass splitting within the isospin 
multiplet, the partial width into the specific charge channel 
should be multiplied by the flavor multiplicity factor ;F (Table 
I Vill i. This 'F factor also incorporates the statistical factor 1/2 
if the final state mesons B and C are identical (as illustrated in 
Figl^. More details of f can be found in the Appendix A of 
Refl64ll. For X(1870) as the mi ^^Di) nh meson, the partial 
decay widths are shown in Table IVinl except the channels of 
S + P mesons. 
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d1 d2 

FIG. 6: qq meson decay diagrams in the ^Pq decay model. 
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TABLE VII: Flavor weight factors. Here, I;;) = {\nn} - \ss))/ V2 and 
I?/') = {\nh) + \ss))/ yl have been taken for simplicity. 



Furthermore, the flux tube model was adopted to check the 
possibility of X( 1870) as a hybrid meson. The partial widths 
are also listed in Table IVIlTl for the comparison. Details of the 
flux model are collected in the Appendix B. 
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TABLE Vni: Predictions for the decay width of X(I870) and compared with results from Ref. 
mode is forbidden. In Ref. 15311 , the mass of hybrid states is taken as 1.8 GeV. 



. A symbol "x" indicates that a decay 



r 



In Ref. 1I23I1 , the strength y is regarded as a free momentum independent constant, and the scale parameter y6 was taken as 



8 



0.4GeV for all mesons. Differently, the momentum depen- 
dent creation vertex has been considered here. The effective's 
(Table m are adopted in our calculations as that in Ref. ll65||. 

Two groups of the partial widths predicted by the Ref. ll53ll 
are quoted in the Table IVIIII The left column was given by 
the flux tube decay model of Isgur, Kokoski, and Paton (IKP) 
with the "standard parameters" [66]. The right column was by 
the flux tube decay model of Swanson-Szczepaniak. 

As shown in Table I Villi our results do not support X(1870) 
as the 772(2'D2) state since its observed decay width is much 
smaller than the the theoretical estimate. Furthermore, the 
BESIII Collaboration claimed that X(1870) primarily decay 
via flo(980)7r 1.17,1 . However, a2(1320)7r is the largest decay 
channel in our numerical results and in Ref. 8231 . If the partial 
width of flo(980);7r channel is amount to a2(1320)7r, the width 
of X( 1870) will be much larger than the observed value. 



ac(??)^ ^« 




X(1870) 



FIG. 7: The diagram for the " <S + P" channels through a virtual 
intermediate I^Pq qq meson. 



In the following, we shall evaluate the partial widths of 
"S + P" channels which have not been given in the tabl dVIIIl 
As illustrated in the Fig|7] we assume X(1870) decay into 
ao(980);r via a virtual intermediate I^Pq qq meson . We 
also notice that 7/(1295) dominantly decay into the rjnn. The 
three-body decay can occur via three intermediate processes: 
7/(1295) ^ 7/(r/flo(980);r/7/(;r;r)s-u«>.. ^ ^^{11 7/(1295) is 
proposed to be the first radical excited state of 7/(550). With 
the ratio r(ao(980)7r)/r(?/7r7r) = 0.65 + 0.10 and r(77(1295)) = 
55 ± 5MeV, the partial width of 7/(1295) to flo(980);r is esti- 
mated no more than 45MeV. The uncertainty of the coupling 
vertex of l^Poiqq) — > ao(980) (see in Fig|7]i is assumed to 



be canceled in the ratio of 



r(X(1870)^ao(980);r) 



r(„(i295)^«o(980w ■ The assumption 
here seems a little rough. However, we just need to evaluate 
magnitudes of these decay channel. 

In our results, the total decay widths of r(X(1870) -^ S + 
P) is no more than 12.6MeV and r(X(1870) ^ flo(980)7r) < 
3.8MeV. The small decay width of r(X(1870) -^ flo(980)7r) 
also indicates that the X(1870) can't be interpreted as the 2'D2 
qq state. 

For a hybrid meson, X(1870) seems most possible to be the 
?//f(0^2"^) state because the total widths exclude the channels 
of S + P are much narrow in our work and in Ref. [1531 . It is 
consistent with the width of X( 1870). 

The assignment for X( 1 870) as the /h(0^ 1 ^^) hybrid seems 
impossible since the theoretical width is rather broad. But 
the widths given by the developed flux tube decay model of 
SS are much small (Table IVlIIl i for the /h(0+1++) hybridllsl. 
So the possibility of X(1870) as a /h(0+1++) hybrid can not 
be excluded. We ague that detecting the decay channel of 



fli(1260)7r can pin down this state. Because this channel for 
the 7///(0^2"^) state is forbidden in the flux tube decay model 
of IKP and very small in the SS flux tube decay model. But it 
is large enough to be observed for the /h(0^1^^) state in two 
kinds of flux tube decay model (see in Tabl jVIlIb . 

As shown in Table IVIIll X(1870) is impossible to be the 
T/zf (O'^O"^) hybrid state. The predicted width in both our work 
and in Ref. |53] are broader In addition, rjn is a visible channel 
for both flo(1450) and 02(1320). A week signal was found in 
the region of 1200~1400MeV in the analysis of tjtt- (Fig. 2(b) 
of Ref.| 17]), which contradicts the large flo(1450)7r channel of 
the riniO^O^^) state. This also implicates X(1870) disfavor a 
0^0"^ hybrid meson. 

In our work and the Ref. ll23ll . the channels of fl2(1320)7r 
and /2(1275)77 are large for the rjiii^Di) state, consistent with 
the observations of the 7/2(1870). However, the partial widths 
of K*K, pp and ww in our work are much narrower than the 
expectations from Ref. f23^. With the mass of 1842 ± 8MeV, 
the decay width of 7/2(1870) is predicted about lOOMeV which 
is much lower than the experiments. 

If X( 1870) and 7/2(1870) are hybrid and quarkonium states 
respectively with the 0^2"^ quantum numbers, they should 
mix with each other Then the interference enhancement will 
enlarge the width of 7/2(1870). The broad decay width of 
7/2(1870) could be explained naturally. On the other hand, 
7/2(1870) has been observed in the channel of flo(980)7r. How- 
ever, this channel seems much small if X( 1870) is a pure 2'D2 
nn meson. The mixing effect will also enlarge this partial 
width. Here, we don't plan to discuss X(1870) and 7/2(1870) 
further in the mixing scheme for the complex mechanism. 



IV. CONCLUSIONS AND DISCUSSIONS 

In this work, we analyze the masses and the two-body 
strong decay of 7/2(1870) and X(1870). The former listed by 
PDG is still in dispute. The latter was observed by the BesIII 
Collaboration recently. Its mass is compatible with 7/2(1870) 
but the decay width is much smaller 

Firstly, the average masses, M„i, are calculated in the the 
Godfrey-Isgur potential model. The presented M„i consist 
with the experimental data well. According to the results, 
X(1870) could be the /;(2Vi), ^\{\^D^) and 7/2(21^2) states. 
In RTs framework, the possible assignments for X(1870) are 
the 7/(3'5o) and f[{l^P\) states. 

Based on the SU(3) flavor symmetry, we extract the 
model independent coupling constants for the processes of 
f{nh) I f {ss) -^ S + P. Then the following conclusions are 
reached: aQji, crrj and /qt/ are the dominant channels for fijifi); 
kK and /o7/ are the only decay channels for f'(ss). For the 
strong decay properties, both 7/2(1870) and X(1870) favor a 
7/2(2'D2) meson assignment if they are a nh meson. 

As the 7/2(2*^2) quarkonium, the two-body strong decay 
of X( 1870) is computed in the ^Pq model. The decay width 
looks much larger than the observation of X( 1870). The broad 
resonance, 7/2(1870), can be a natural candidate for the 2'D2 
nfi meson. There, we fix the creation strength, y, in two kinds 
of processes: ®.nh -^ nh + nh; ®. nh ^ ns + ns and ss — > 



ns + ns. The functions of creation vertex are determined as 
yip) = 3.0 + 25exp(-V) and y{p) = 1.8 + Aexp(-lOp^) 
respectively. 

For the unclear nature of the S mesons, we just evaluate the 
magnitude of the partial widths of ">S + P" channels by the 
ratio, r(^(i295)^fl°(980)^) . under a rather crude assumption that 
^(1295)/X(1870) -^ flo(980)7r through a virtual intermediate 
l^Po qq meson (see Fig|7]). The uncertainties of the coupling 
vertex for l^Po(qq) -^ ao(980) are assumed to be canceled in 
the ratio. The total widths of "S + P" are evaluated no more 
than 7.3MeV. Since X(1870) primai'ily decay via flo(980)7r, 
it also indicated that the X(1870) can't be interpreted as the 
2^D2 nh state. 

If X(1870) is a hybrid states, the 7^'' quantum numbers 
must be 0^0"^, 0^1^^ and 0^2"^. Here, we employ the flux 
model to study the hybrid possibilities. Our results agree well 
with most of predictions given by Ref. |53]. X(1870) seems 
most possible to be the 77//(0^2"^) state, because the predicted 
widths are quite narrow, which is consistent with the width 
of X(1870). But we can'texclude the possibility of 0+1++. A 
precise measurement of a i ( 1 260)7r and 02(1 320)7r is suggested 
to pin down the underlying structure of X(1870). 

In addition, the mixture shall occur between 772(1870) and 
X(1870) when they have the same j'''^ of 0^2"^. If the mixing 
occurs, the interference enhancement should enlarge the width 
of 772(1870), which can explain the broader decay width of 
7/2(1870) naturaUy. 

In summary, we propose that there are two diff'erent 
isoscalar mesons around 1.8GeV. The broad one, 7/2(1870), 
can be interpreted as the 7/2(2*^2) quarkonium. The narrow 
one, X(1870), looks more like the 7/2(0^2"^) hybrid meson. 
The possibiHty of X( 1870) as a ///(O+T^) hybrid still exists. 
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Appendix A: The expressions of amplitudes 

We have omitted a exponential factor in following decay 
amplitudes Mls for compactness. 
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The above amplitudes, Mls , can be reduced further in the 
approximation of mi - mz - m and Pa - Pb - Pc - P - 
The reduced Mls are consistent with these given by Ref. ll64ll 
except for an unimportant factor, -2^/^, since this factor can 
be absorbed into the coefficient y. 



Appendix B: Hybrid decay In the flux tube model 



The flux tube model was motivated by the strong coupling 
expansion of the lattice QCD. In this model, decay occurs 
when the flux-tube breaks at any point along its length, with a 
qq pair production in a relative 7^'' = 0^^ state. It is similar 
to the ^Po decay model but with an essential difference. The 
flux tube model extend the nonrelativistic constituent quark 
model to include gluonic degrees of freedom in a very simple 
and intuitive way, where the gluonic field is regarded as tubes 
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of color flux. Then it can be extended to the hybrid research. 
When the hybrid mesons are assumed to be narrow, and the 
threshold effects aren't taken into account, the partial decay 
width Fls (H — > BC) is given by the flux model as ll67ll 



where Ma, Mb, Mc are the "mock-meson" masses of A, B, 
C 1 65]. When a hybrid meson decay into f-wave and pseu- 
doscalar mesons, the partial wave amplitude Ml{H -^ BC) 
(with S = 5 b) is given as the following form 



Tls(H^BC) 



MbMc 



{2JA + l)n Ma 



\MLs(H^BCr (Bl) 



MdH ^ BC) = {(pEMcpAMi 



9V32 «W;r\l 
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„3/2+,5 

I-'a 



•/fe/(2^2))2 V3r(3/2 + 5) p 



MdH -^ BC) 



(B2) 



r 



B 



0++ 

2++ 



0+0-+(H) 



■V2Ms/3 



+Mo/3 



0+l++(H) 



-V2Mp,/V3 
-Mp,/V2 
-Mp J ^/30 



0+2-+(H) 



+Md/3 

-VSMs/VTS 
-V7Mo/3 



TABLE IX: Partial wave amplitudes Ml(H -> BC) for an initial 
hybrid H decaying into a P-wave and pseudoscalar mesons. 



The flavor matrix element {(pB<Pc\<pA<po} have been discussed 
before. Ml(H -^ BC) are listed in Table |IX] for the states of 
miO^O'^), MOn^^) and mio^i-^). 

Here the Ms, Md, Mp. and Mr are defined as Ms - 
-(3/io-|i +4/12), VHd = (|i +5/^2), W(f, =-/(2|o + 3/ii-|2), 
Mp, = -i{go + §2), Mp, = -/(lOgo + 9/Ji + I2) and Mf = 
-3i(g2 + hj,). The analytical expressions of §,■ and /i, are given 
as 



J 



= 2 



3+d 



M"m 



(M + my 



,^^7 7,7n "+*+3 ^ n + 5 + 3 , n + 6 + 3 , , M 

-{ipl+P^—Yi )yFd ,«+l,-( 



M + m 



2/32+^2 



^^P" 



(B3) 
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M + m 



)"i2/3l+^^) 



,ts+t n + 6 + 4 n + 6 + 4 

~n — - — )iFA — - — ,« + 1, 



M 



M + m 2I3\+P' 



■]P" 



r 



(B4) 



where iFi[-- ■] are the confluent hypergeometric functions. 
Here we don't take account of the decay channels of H ^ 
2S + IS because they are forbidden by the conservation laws, 
or the "spin selection rule", or the phase space, e.g., the decay 
channel of 7r( 1 300) + tt is forbidden for the /h(0^ 1 ^^) state by 
the "spin selection rule". In this work, we choose to follow 

the Refs. ||67D and take the combination (ac/9 V3)5A°) '-^ 



as 0.64 which was fixed by the conventional mesons II65I1 . 
M ^ m ^ mu,d = 330MeV, M'g=^ = M'=^ = 1250MeV, 
M'=° = 720MeV, M'=^ = 850MeV,y6^ =0.27GeV, S =0.62, 
b =0.18GeV2 and k =0.9. Final states containing n have 
^ =0.36GeV, otherwise j8 =0.40GeV. 



[1] J. Beringer et al. (Particle Data Group), Phys. Rev. D 86, 

010001 (2012). 
[2] C.A. Meyer and Y. Van Haarlem, Phys. Rev. C 82, 025208 

(2010). 
[3] B. Aubert et al. (BABAR Collaboration), Phys. Rev. D 76, 

012008 (2007). 
[4] J. P Lees et al. (BABAR Collaboration), Phys. Rev. D 86, 

012008 (2012). 
[5] M. Ablikim et al. (BBS Collaboration), Phys. Rev. Lett. 100, 



102003 (2008). 
[6] C. P Shen et al. (Belle Collaboration), Phys. Rev. D 80, 

031 lOl(R) (2009). 
[7] Gui-Jun Ding and Mu-lin Yan, Phys. Lett. B 650, 390 (2007). 
[8] Gui-Jun Ding and Mu-lin Yan, Phys. Lett. B 657, 49 (2007). 
[9] Zhi-Gang Wang, Nucl. Phys. A791, 106 (2007). 
[10] Hua-Xing Chen, Xiang Liu, A. Hosaka and Shi-Lin Zhu, Phys. 

Rev. D 78, 034012 (2008). 
[11] Shi-Lin Zhu, Int.J.Mod.Phys.E 17, 283, (2008). 



12 



[12] K. Karch et al. (Crystall Ball Collaboration), Phys. Lett. B 249, 

583 (1990). 
[13] M.Feindt,DESY-90- 128 (1990). 
[14] J.Z. Bai et al. (BBS Collaboration), Phys. Lett. B 446, 356 

(1999). 
[15] J. Adomeit et al. (Crystal Barrel Collaboration), Z. Phys. C 71, 

227 (1996). 
[16] A.V. Anisovich et al., Phys. Lett. B 477, 19 (2000). 
[17] M. Ablikim et al. (BBS Collaboration), Phys. Rev. Lett. 107, 

182001 (2011). 
[18] M.S. Chanowitz and S.R. Sharpe, Nucl. Phys. B222, 211 

(1983). 
[19] T. Barnes, P. Close and F. de Viron, Nucl. Phys. B224, 241 

(1983). 
[20] N. Isgur and J.B. Paton, Phys. Rev. D 31, 2910 (1985). 
[21] T. Barnes, P.B. Close and B.S. Swanson, Phys. Rev. D 52, 5242 

(1995). 
[22] S. Ishida, H. Sawazaki, M. Oda and K. Yamada, Phys. Rev. D 

47, 179 (1993). 
[23] De-Min Li and Bn Wang, Bur Phys. J. C 63, 297 (2009). 
[24] S. Godfrey and N. Isgur, Phys. Rev. D 32, 189 (1985). 
[25] E. Bichten, K. Gottfried, T. Kinoshita, J.B. Kogut, K.D. Lane 

and T.M. Yan, Phys. Rev. Lett. 34, 369 (1975). 
[26] B.B. Salpeter and H.A.Bethe, Phys. Rev. 84, 1232 (1951). 
[27] B.B. Salpeter, Phys. Rev. 87, 328 (1952). 
[28] D.S. Hwang, C.S. Kim and W. Namgung, Z. Phys. C 69, 107 

(1995). 
[29] S.R Radford, W.W. Repko and M.J. Saelim, Phys. Rev. D 80, 

034012 (2009). 
[30] S.R Radford and W.W. Repko, Phys. Rev. D 75, 074031 (2007). 
[31] W. Roberts and B. Silvestre-Brac, Phys. Rev. D 57, 1694 

(1998). 
[32] T. Barnes, S. Godfrey and B.S. Swanson, Phys. Rev. D 72, 

054026 (2005). 
[33] A.M. Badahan and B.L.G. Bakker, Phys. Rev. D 64, 114010 

(2001). 
[34] A.M. Badalian, B.L.G. Bakker and Yu.A. Simonov, Phys. Rev. 

D 66, 034026 (2002). 
[35] P Gonzalez, Phys. Rev. D 80, 054010 (2009). 
[36] S. Godfrey and R. Kokoski, Phys. Rev. D 43, 1679 (1991). 
[37] RB. Close and B.S. Swanson, Phys. Rev. D 72, 094004 (2005). 
[38] De-Min Li, Peng-Rei Ji and Bing Ma, Bur. Phys. J. C 71, 1582 

(2011). 
[39] D. Bbert, R.N. Baustov and V.O. Galkin, Phys. Rev. D 79, 

114029(2009). 
[40] R. Giachetti and B. Sorace, Phys. Rev. D 87, 034021 (2013). 



[41] P. Masjuan, B.R. Arriola and W. Broniowski, Phys. Rev. D 85, 

094006 (2012). 
[42] L. Maiani, P. Piccinini, A.D. Polosa and V. Riquer, Phys. Rev. 

Lett. 93, 212002 (2004). 
[43] J.D. Weinstein and N. Isgur, Phys. Rev. Lett. 48 , 659 (1982); 

Phys. Rev D 27, 588 (1983); Phys. Rev. D 41, 2236 (1990). 
[44] R.L. Jaffe, Phys. Rev. D 15, 281 (1977); R.L. Jaffe and RB. 

Low, Phys. Rev .D 19, 2105 (1979); M. Alford and R.L. Jaffe, 

Nucl. Phys. B578, 367 (2000). 
[45] Zhi-Gang Luo, Xiao-Lin Chen and Xiang Liu, Phys. Rev . D 

79, 074020 (2009). 
[46] Jie-Sheng Yu, Zhi-Peng Sun, Xiang Liu and Qiang Zhao, Phys. 

Rev. D 83, 114007(2011). 
[47] Jia-Peng Liu, Gui-Jun Ding and Mu-Lin Yan, Phys. Rev. D 82, 

074026 (2010). 

[48] R.L. Jaffe, Phys. Rept. 409, 1 (2005)^ 

[49] A. Selem and R Wilczek, hep-ph/06" 02T28] 

[50] B. Santopinto and G. Galata, Phys. Rev:"(r75, 045206 (2007). 

[51] J.J. de Swart, Rev. Mod. Phys. 35. 916 (1963). 

[52] C. Amsler and RB. Close, Phys. RevD 53, 295 (1996). 

[53] PR. Page, B.S. Swanson and A.P. Szczepaniak, Phys. Rev. D 

59, 034016 (1999). 
[54] L. Micu, Nucl. Phys. BIO, 521 (1969). 
[55] A. Le Yaouanc, L. Olivier, O. Pene, and J.C. Raynal, Phys. Rev. 

D 8, 2223 (1973); A. Le Yaouanc, L. Oliver, O. Pene, and J.C. 

Raynal, Hadron Transition in the Quark Mo(ie/(Gordon and 

Breach, New York, 1988), p. 311. 
[56] Zao-Chen Ye, Xiao Wang, Xiang Liu and Qiang Zhao, Phys. 

Rev. D 86, 054025 (2012). 
[57] C. Hayne and N. Isgur, Phys. Rev D 25, 1944 (1982). 
[58] M. Jacob and G.C. Wick, Ann. Phys. (N.Y) 7, 404 (1959). 
[59] P Geiger and B.S. Swanson, Phys. Rev. D 50, 6855 (1994). 
[60] H.G. Blundel, arXiv:hep-ph/9608473 
[61] R. Bonnaz and B. Silvestre-Brac, Pew-Body Syst. 27, 163 

(1999). 
[62] B.S. Ackleh, T. Barnes and B.S. Swanson, Phys. Rev. D 54, 

6811 (1996). 
[63] J. Segovia, D.R. Bntem and B. Bemandez, Phys. Lett. B 715, 

322 (2012). 
[64] T. Barnes, RB. Close, PR. Page and B.S. Swanson, Phys. Rev. 

D 55, 4157 (1997). 
[65] R. Kokoski and N. Isgur, Phys. Rev. D 35, 907 (1987). 
[66] N. Isgur, R. Kokoski, and J. Paton, Phys. Rev. Lett. 54, 869 

(1985). 
[67] RB. Close, PR. Page, Nucl. Phys. B443, 223, (1995). 



